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Abstract 

We consider the evolution of a flat Fricdniann-Robcrstson- Walker Universe in a 
higher derivative theory, including a terms to the Einstcin-Hilbert action in the 
presence of a variable gravitational and cosmological constants. We study here the 
evolution of the gravitational and cosmological constants in the presence of radiation 
and matter domination era of the universe. We present here new cosmological solutions 
which are physically interesting for model building. 
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I. INTRODUCTION : 

In recent years there has been a growing interest to study cosmological models with vari- 
able cosmological constant [1] and gravitational constant [2]. The interest stems from the 
observational analysis of type la supernova with red shift parameter Z < 1 [3] . Supernovae 
study provides evidence that we may be living in a low mass density universe of the order of 
n ~ 0.3 [4]. However, inflationary universe model which is an essential ingredient in mod- 
ern cosmology predicts, density parameter Q 1. Thus a major part of the matter content 
in the universe remains unobserved, which we call dark matter. It is also understood from 
the observational analysis of type la supernovae that the present universe is accelerating. 
This leads to the assumption that there exists a small positive cosmological constant (A) 
even at the present epoch. It is also understood in modern Cosmology that our universe 
might have came through a phase of rapid expansion or inflation in the past. Such a rapid 
expansion or inflation may be obtained if one considers an Einstein-Hilbert action with a 
cosmological constant. To accommodate and to interpret the above facts several physical 
models are proposed in recent times with variable cosmological and gravitational constant. 
The effect of extra dimensions in the evolution of a homogeneous isotropic universe and a 
first quantitative study on the space and/or time dependence of fine structure and New- 
ton's constants were taken up by Forgacs and Horvath [5]. It was shown that the time 
and space variation of fundamental constants is a natural consequence of Kaluza-Klein 
theories. Bertolami [6] studied cosmological model with a time dependent cosmological 
term. Later, Ozer and Taha [7] and others [8] obtained models of the universe, which are 
free from cosmological problems considering a varying cosmological constant in the case of 
a closed universe. Recently, a number of literature [9-17] appeared which describe various 
cosmological issues, for example, the age problem, observational constraints on A, struc- 
ture formation and gravitational lenses. Ratra and Peebles [10] studied the dark energy 
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content of the universe considering a varying cosmological constant which may be reahzed 
using a scalar field. Dolgov [11] analysed the cosmological evolution of a free massless vec- 
tor or tensor (not gauge) field minimally coupled to gravity and found that cosmological 
constant exactly cancels out. Sahni and Starobinsky [12] discussed the observational and 
theoretical aspects of a small positive cosmological constant term which may be relevant 
at the present epoch. The origin of such a small A has been discussed in the framework of 
field theoretic techniques and by modelling a dynamical A term by a scalar field. Recently, 
Padmanabhan [13] disscussed several aspects of the cosmological constant both from the 
cosmological and field theoretic perspectives in order to explain the cosmological constant 
problem in the framework of the string theory. Vishwakarma [14] explored different models 
to estimate the density parameter and/or deceleration parameter for the variable A models 
and numerically analysed the observational data in cosmology. It is found that the models 
with dynamical cosmological term (A) are acceptable with observations. Birkel and Sarkar 
[15] studied the nucleosynthesis and obtained a bound on a time varying A. 

In the last few years a number of work reported in cosmology with variable cosmological 
and gravitational constant with perfect fluid as a source in isotropic [16] as well as in 
anisotropic models [17] of the universe in the frame work of Einstein gravity. In this paper 
we study isotropic cosmological models in a higher derivative theory (which is also known 
as a generalized theory) with variable cosmological and gravitational constant to explore 
cosmological solutions which can accommodate the recent observations in the universe. A 
number of literature appeared in the last few decades on the generalized theory of gravity 
with additional terms in R^, Rfj^i.psR'^^^^ , C^npsC^'^'^^ in the Lagrangian. However, in 4 
dimensions only R^ -term is important and the combinations of terms e.g. Gauss-Bonnet 
combinations occurs as Euler number in the Einstein action. The generalized theory of 
gravity has a number of good features. It is renormalizable, asymptotically free theory 
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[18]. However, the unitarity problem of the theory is not yet settled [19]. 

The paper is organized as follows : in sec. II, we present the gravitational action and 
field equations, Cosmological solutions are given in sec. Ill and a brief discussions in sec. 
IV. 



II. GRAVITATIONAL ACTION AND FIELD EQUATIONS : 

We consider a gravitational action with higher order terms in the scalar curvature (R) 
with a variable cosmological constant (A(t)) and variable gravitational constant (G(t)) 
which is given by 



fiR) 



+ Lr, 



— j4 

-g a X 



(1) 



.167rG(t) 

where f{R) is a function of R and its higher power including cosmological constant and 

Lm represents the matter lagrangian. 

Variation of action (1) with respect to g^n yields 



f'{R) R^, - - f{R) g^, + f"{R) {V^V.R - URg^,) + 



f"'{R) {V^RV.R - V'RVaR gt..) = - SttG T^, 



(2) 



where □ = g^u'V^'V'^ and is the covariant differential operator, and prime represents 
the derivative with respect to R, T^^ is the energy momentum tensor for matter determined 
by Lm . 

We consider flat Robertson- Walker spacetime given by the metric 



ds-" = -dt^ + a\t) dr^ + r^{de^ + sin^ed(j) 



2/iJ 



where a(t) is the scale factor of the universe . The scalar curvature in this case is 



(3) 



R = -Q[H + 2H'^ 



(4) 
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where H = ^ is the Hubble parameter. The trace and (0,0) components of eq.(2) are given 
by 

Rf'iR) - 2f{R) + Sf"{R) (^R + 3^i?^ + 3f"\R)R + SttG T = 0, (5) 

f'{R)Roo + h{R)-Sf"{R)-R + 8TrG Too = (6) 

where dot represents derivative with respect to comoving time (t). 

Now we consider f{R) = R + aR^ — 2A{t). The equations (5) and (6) yield 



H'^-6a 



2HH -H'^ + 6HH'^ 
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where p is the cosmic fluid energy density and p represents the pressure. The equation of 
state is 

P = (7 - 1)P (9) 
where 7 is a constant (1 < 7 < 2). The deceleration parameter (g) is releted to as 



dt \H, 

III. COSMOLOGICAL SOLUTIONS : 

We now use equations (7) — (9) to obtain cosmological solutions. The set of equations 
contain six unknown, to solve it we apriori consider time variation of two of them for a 
given equation of state of the early universe. 

Let us consider power law behavior of the universe which is given by 

a{t) = ao (10) 

where oq and n are the constants. We now look for cosmological solutions in two different 
cases: (i) Radiation dominated and (ii) Matter dominated regime of the universe respec- 
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tively in the next sections 



(i) Radiation Dominated Universe : 

In this case we have 7 = | i-e. p = \p and in the absence of particle creation eq.(8) can 
be written as 

^ + AHp = 0, (11) 
A{t) + 8TrpG{t) = 0. (12) 
Eq.(ll) yields the energy density with the scale factor of the universe, which is 

.= 0. (13) 

where po is a constant. Using eqs.(7),(12) and (13), one gets evolution of cosmological and 
gravitatio; S. Perlmutter, Astro. Phys. J. 517, 517 (1999)nal constants, which are given 
by 

A(t) = |l(2n-l) + ^^""("-;^(^"-^\ (14) 



G{t) 



3d, 



4 



Wirpo 



n 36an{2n — 1) 



+ 



(15) 



^2— 4ti ^4— 4n 

The deceleration parameter becomes: q = which is positive for < 1 and 

negative for n > 1. 

We note the following cosmological solutions : 
(i) When n=\ i.e. a ~ \/t, we get 

A = ; g > 

and 

3an 

G = = constant. 

SzTTpo 

Thus it is evident that the presence of higher derivative terms do not effect the behaviour 
of the early universe in radiation dominated universe, in the absence of a cosmological 
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constant and a constant G. It has a particle horizon. The cosmological evolution in this 
case is the same as that obtained for a radiation dominated era in the Einstein gravity, 
(ii) When n = 1 i.e. a ~ t, we get 



and 

G{t) 



3oo 



+ 36q; 



In this case k{t) decreases with time where as G{t) increases. To begin with G is a con- 
stant and it is determined by the coupling constant (a), initial size and matter density. 
For a < 0, the solution is physically relevant for the epoch t > 6-\/|a|. In this case the 
universe has no particle horizon. 



We now discuss the solution permitted by the field equations in the presence of particle 
creation during radiation domination. The eq.(8) with radiation becomes 



dp 
'dt 



+ 4pH=- 



G A 
G^^8^ 



(16) 



The field eqs.(7) and (16) give the energy density and the cosmological constant which are 

3n 



167rG 



1 36a(2n - 1) 
i2 + ^4 



(17) 



and 



3n , , 54an(2n- l)(n- 1) 
A = ^(2n-1) + 



(18) 



We now discuss special cases : 
(i) When n = ^, we get 



A = 0, 

3 1 
" 32TrG¥ 
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Here the solution is similar to the usual matter dominated universe in the Einstein gravity, 
which we get even if a 7^ in the absence of cosmological constant, 
(ii) When n = 1, we get 

, 3 



P 



1 36a 



For n = 1, the cosmological constant is independent of the higher derivative term, however, 
the energy density is determined by a coupling parameter a of the higher derivative theory. 
For a < 0, physically relevant solutions are obtained for t > 6-v/|a[. The above solutions 
are permitted with a gravitational constant which remains invariant. 



(ii) Matter Dominated Universe: 

In this case we consider p = {) i.e. 7 = 1, and in the absence of particle creation eq.(8) can 

be written as two decoupled equations, given by 

^ + 3Hp = 0, (19) 
A{t) + 87rpG{t) = 0. (20) 

Eq.(19) yields the energy density 

(21) 

where po is a constant. Using eq.(7),(20) and (21), the cosmological constant and gravita- 
tional constant are obtained, which are 

n, „, 18Q!n(2n — l)(3n — 4) 
Ait) = -(3n - 2) + ^ ^ 1, (22) 



and 



G{t) 



AlTpQ 



n ^ 36an(2n — 1) 



f2-3n ' t4-3n 



(23) 
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We note the following cosmological solutions : 



(i) When n 



I.e. a ~ 



we get 



and 



G{t) 



q>0 



In this case one gets solution with negative cosmological constant , however, the Gravi- 
tational constant remains positive and decreasing with time. Thus a comparatively less 

2 

expansion in the matter dominated regime is obtained in this case ( since a(t) ~ ts, in 
Einstein gravity). 



(ii) When n = | i.e. a ~ , we get 



q>0 



and 



G{t) 



47rpo 



2 8a' 

3 ^ 



We obtain Matter dominated solution with A < if a > and A > if a < 0, however 
G attains a constant value at a later epoch, i.e. at t ^ 00. It has a particle horizon. 
This solution is new and interesting in model building for describing late universe. The 
cosmological constant becomes small depending upon time and the coupling parameter at 
late epoch. Here one recovers the matter dominated solution of Einstein gravity for a = 0,( 
as A = and G = Go). 

(iii) when n = 3, i.e. a ^ , which is a rapid power law expansion. In this epoch, we get 



m = 



q<0 



and 



G{t) 



80a + -r 
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Thus in the matter dominated case an accelerating universe results in the presence of 
a decaying cosmological constant, but with an increasing gravitational constant G{t) for 
a > 0. For a < 0, at t — 0, one gets solution with G{t) negative which later attains 



positive values for t > y^|60a| and thereafter G{t) grows. The solution has no particle 
horizon. 

In the presence of particle creation and matter domination era, eq.(8) becomes 

Prom eqs.(7) and (24), the energy density and the cosmological constant are obtained, 
which are 

n ri 36Q;f2n — 1)1 ^^^^ 



AttG 



1 ^ 36a(2n - 1) 



t2 t4 



18an(2n- l)(3n-4) 
A=-{3n-2) + ^ ^ L (26) 

We now discuss special cases : 

(i) when n = ^ i.e. a ~ ^/i, the cosmological parameter and the energy density varies as, 

_ 1 

(ii) When n = | i.e. a ^ ta , we get 

8a 

^ = T, 



6ttG 



1 12a 

t2 



Thus the usual matter dominated solution is recovered with A < here. 

(iii) when n = 3, i.e. a ~ ts, which is a rapid power law expansion. In this epoch, the 

cosmological constant and the energy density become 

A = 



10 



1 



1 60a 



i2 ^4 



For n = ^ and n = |; the cosmological constant is independent of the higher derivative 
term, however, in the later case energy density also depends on the coupling parameter 
associated with the higher derivative term. In the later case it has no particle horizon. 
But in the case of n = |, both the cosmological constant and energy density is determined 
by a. It may be pointed out here that the gravitational constant G becomes a constant in 
this case. The solution obtained here describes a universe which is still accelerating in the 
matter dominated phase. 
IV. Discussions : 

In this paper we investigate isotropic cosmologies with varying G and A in a i?^ theory. 
We studied both radiation dominated and matter dominated phases of evolution of the 
universe with varying G and A and obtained some new cosmological solutions, which may 
be important to describe a late universe. We also discussed relevance of some of the so- 
lutions as special cases. We obtain here some solutions where G is increasing similar to 
that found by Abdel Rahaman [2] and Arbab [20] in Einstein gravity. The cosmological 
constant A, however, decreases with time in those cases. In the frame work of theory 
we note an interesting solution during matter dominated regime. The gravitational con- 
stant G decreases initially and ultimately it attains a constant value at a later stage of 

2 

evolution of the universe. The universe evolves as a ~ which is similar to that one gets 
in the Einstein gravity during matter domination but we note a realistic variation for the 
parameter G. In this scenario G may be large in the early universe, however, at the present 
epoch it becomes a constant. The above solutions are permitted either with (i) A > if 
a < or (ii) with A < if a > . We also note that in the matter dominated era an 
interesting new solution exists in theory which admits an accelerating universe (sup- 
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ported by supernova observation), with a decaying cosmological constant and a constant 
G. It may be pointed out here that an inflationary solution (a(i) ~ e^°*) are permitted in 
-theory if one considers a non zero A even with p = p = 0. Power law inflation is also 
admitted in the matter dominated universe (MDU) or in the radiation dominated universe 
(RDU) with varying A and G. In MDU, power law expansion of the universe is permissible 
if one allows G to be increased at a later stage. However, when A — > 0, the solution is 
no longer valid and the usual matter dominated universe solution may be recovered with 
a constant G. The stability analysis of the cosmological solutions discussed here will be 
taken up elsewhere. 

Ackknowledgements : 

The authors would like to thank lUCAA Reference Centre at North Bengal University 
for providing facilities during the work. BCP would like to thank University Grants Com- 
mission, New Delhi for awarding Minor Research Project (Project No. 30-16/2004(SR)) 
for completing the work. 



12 



References 

[1] S. L. Adler, Rev.Mod.Phys. 54, 729 (1982) and references therein. 

[2] A-M M. Abdel-Rahaman , Gen. Rel. Grav. 22, 655 (1990); A. Beesham, Nuovo 
Cimento B 96, 17 (1986); A. Abdulssattar and R. G. Vishwakarma, Class Quantum 
Grav 14, 945 (1997); A. I. Arbab, Gen. Rel. Grav. 29 391 (1998); Astrophys J. 565, 
483 (1997); A. G. Riess et al, Astron. J. 116, 1009 (1998). 

[3] S. Permutter, et al. Nuture 51, 391 (1998) ^Astrophys J.565, 483 (1997); A. G. Riess, 
et al. Astron. J. 116, 1009 (1998); S. Perlmutter, Astro. Phys. J. 517, 517 (1999). 

[4] A. G. Riess, Nuovo. Cim. B 93, 36 (1986). 

[5] R Forgacs and Z. Horvath, Gen. Rel. Grav. 10, 931 (1979) ; P. Forgacs and Z. Horvath, 
Gen. Rel. Grav. 11, 205 (1979). 

[6] O. Bertolami, Astron. J. 116, 1099 (1998) . 

[7] M. Ozer and M. O. Taha, Phys. Lett. B 171, 363 (1986); Nucl. Phys. B 287, 776 
(1987). 

[8] S. Weinberg, Rev. Mod. Phys., 61, 1 (1989) ; A-M M. Abdel-Rahaman, Phys. Rev. 
D 45, 3497 (1992); W. Chen and Y. S. Wu, Phys. Rev. D41, 695 (1990); Erratum, 
Phys. Rev. D 45, 4728 (1992). 

[9] W. L. Freedman, Phys. Scripta T85, 37 (2000); A. Balbi et al , Astro. Phys. J. LI, 
545(2000). 

[10] P. J. E. Peebles and B. Ratra, Ap. J. 325, L17 (1988); Phys. Rev. D 37, 3406 (1988) 
[11] A. D. Dolgov, Phys. Rev. D 55, 5881 (1997). 



13 



[12] V. Sahni and A. Starobisky, Int. J. Mod. Phys. D 9, n the presence of a negative 
cosmological constan373 (2000); astro-ph/9904398; V. Sahni, Phys. Rev. D 42, 453 
(1990). 

[13] T. Padmanabhan, Phys. Rep. 235, 380 (2003). 

[14] R. G. Vishwakarma, Gen. Rel. Grav. 197, 33 (2001); Class. Quantum Grav. 19 , 4747 
(2002); Class. Quantum Grav. 18 , 1159 (2001); R. G. Vishwakarma and Abdussattar, 
Phys. Rev. D 60, 063507 (1999). 

[15] M. Birkel and S. Sarkar, Astroparticle Physics. 6, 197 (1997). 

[16] D. Kahigas, R S. Wegson, C. W. F. Everitt, Gen Rel.Grav. 24, 351 (1992). 

[17] A. Beesham, Gen.Rel.Grav. 26, 159 (1994); D. Kahigas, RS. Wegson, C. W. F.Everitt, 
Gen. Rel. Grav. 27, 3645 (1995); A. Rradhan and A. Kumar, IntJ. Mod. Phys. DIO, 
291 (2001). 

[18] K. S. SteUe, Phys. Rev. D33, 2756 (1986). 

[19] I. Antoniadies and E. T. Tombouhs, Phys. Rev. D 33, 2756 (1986). 
[20] A. I. Arbab, Gravitation and Cosmology, 8 , 227 (2002). 



14 



